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Abstract 

We present the two-loop corrected operator matrix elements calculated in 
iV-dimensional regularization up to the finite terms which survive in the limit 
e = N — 4 — > 0. The anomalous dimensions of the local operators have been 
previously extracted from the pole terms and determine the scale evolution of 
the deep inelastic structure functions measured in unpolarized lepton hadron 
scattering. The finite e-independent terms in the two-loop expressions are 
needed to renormalize the local operators up to third order in the strong 
coupling constant a s . Further the unrenormalized expressions for the two- 
loop corrected operator matrix elements can be inserted into specific one 
loop graphs to obtain a part of the third order contributions to these matrix 
elements. This work is a first step in obtaining the anomalous dimensions 
up to third order so that a complete next-to-next-to-leading order (NNLO) 
analysis can be carried out for deep inelastic electroproduction. 

PACS numbers : 11.15.Bt, 12.38.-t, 12.38.Bx 
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I. INTRODUCTION 



One of the most important successes of the theory of perturbative Quantum Chromody- 
namics (QCD) is the prediction of the scale evolution of the structure functions Fi(x,Q 2 ) 
measured in deep inelastic lepton-hadron scattering (for reviews see Jl|). From these struc- 
ture functions one can infer the parton densities, which serve as input for many other deep- 
inelastic (hard) processes, and lead to a wealth of predictions for cross sections (for a review 
see [0). Many of these predictions have been tested by now and are in very good agreement 
with the experimental data. The scale evolution of the structure functions is determined by 
the anomalous dimensions of local composite operators 0™(x), i = q,g, where n denotes the 
spin, which show up in the expansion of the product of two local electroweak currents J^x) 
and J u {y) around the lightcone (x — y) 2 = 0. The structure functions are proportional to the 
Fourier transform into momentum space of the product of these two currents sandwiched 
between hadronic states. In this way the structure functions can be written as a product of 
operator matrix elements (OME's) and coefficient functions. The former describe the long 
distance (low momenta) properties of QCD whereas the latter account for the short distance 
(large momenta) behaviour of the specific quantities, like e.g. the structure functions, under 
consideration. Both the anomalous dimensions and the coefficient functions are calculable 
order by order in perturbation theory so that they can be expressed in a series expansion 
in the strong coupling constant a s . The finite anomalous dimensions are responsible for the 
scale dependence of the structure functions and the parton densities. On the other hand the 
OME's themselves are of a non-perturbative nature, and, apart from some attempts using 
lattice gauge field theory, cannot yet be determined from first principles. Hence they have 
to be determined from models and fits to experimental data. In the past one has put a 
lot of effort into calculating the lowest order coefficients in the perturbation series for the 
anomalous dimensions and the coefficient functions. In || and [|J the order a s contributions 
to the anomalous dimensions were obtained by evaluating the one-loop OME's which follow 
from inserting the local operators Of(x) between quark and gluon states (for an alternative 
method see ||). The second order contributions, which involved the evaluation of two- loop 
OME's, have been computed by various groups |J- ||. As far as the coefficient functions are 
concerned they have been evaluated for many processes up to next-to-leading order (NLO) 
(for a recent review see M). Combining them with the above NLO anomalous dimensions 
one can make a complete NLO analysis of many quantities like cross sections for Drell-Yan 
production or deep inelastic structure functions. 

The advent of HERA opened up a new era of deep inelastic lepton-hadron experiments 
with much higher statistics than was previously available. Furthermore the values of x and 
Q 2 , on which the structure functions depend, could be extended beyond those accessible by 
earlier fixed target experiments. This enables QCD to be tested with even higher degree of 
precision so that effects beyond NLO can be studied. In particular both the small and the 
large x- regions have attracted much attention in the literature (for a review see e.g. [[Hfl). 
Hence it will be necessary to extend the existing expressions for the anomalous dimensions 
and the coefficient functions beyond NLO. The first step to obtain the next-to-next-to- 
leading order (NNLO) corrections has been made in fly]] where the order a 2 corrections to 
the coefficient functions for deep inelastic lepton-hadron scattering have been computed. The 



same program has also been completed for the Drell-Yan process in [12|. For a consistent 
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NNLO analysis we also need the third order contributions to the anomalous dimensions. 
Unfortunately the latter are not known yet except for those corresponding to the operators 
Of(x) with spin n = 2,4,6,8 and 10. They have been computed in J0| using completely 
different methods than those used in |J- |J. Using the above second order coefficient 
functions and the third order anomalous dimensions it is possible to make a NNLO analysis 
of the structure functions F 2 (x, Q 2 ) and F 3 (x, Q 2 ) as long as one limits oneself to the large 
x (x > 0.01) and small Q 2 (see [Q). If one wants to study the behaviour of the structure 
functions at small x or at large Q 2 knowledge of the anomalous dimensions for any spin is 
indispensable. Since the full calculation of the latter quantities is a tremendous enterprise 
it has to be carried out step by step. The first step is to compute the two-loop OME's up 
to finite terms which is an extension of the work done in ||- ||. These finite terms are 
needed to carry out the renormalization of the three-loop graphs since they determine the 
single pole terms in A-dimensional regularization from which one has to extract the third 
order contributions to the anomalous dimensions. Moreover by inserting the unrenormalized 
two-loop corrected OME's in one-loop graphs one gets parts of the expressions for the 
unrenormalized three-loop graphs. This method has been used previously to get parts of 
the two-loop expressions by inserting the one-loop corrected OME's into one-loop graphs. 

In section 2 we will give an outline of the calculation which will be carried out in the 
Feynman gauge. Although this gauge entails some complications due to the mixing of 
physical and unphysical operators it is the only one in which it is feasible to perform a 
calculation of Feynman graphs for the OME's beyond two-loop order. All other gauges, like 
the axial gauge used in ||, lead to even more complications, which we want to avoid. The 
long expressions obtained for the full OME's are presented in Appendix A while Appendix B 
contains results for non-physical OME's, which are needed to carry out the renormalization 
programme. 



II. THE CALCULATION OF THE TWO-LOOP OPERATOR MATRIX 

ELEMENTS 



In this section we will give an outline of the calculation of the OME's up to two-loop 
order. The operators, which appear in unpolarized lepton-hadron scattering, can be divided 
into singlet and non-singlet parts with respect to the flavour group. In leading twist (namely 
twist two) the non-singlet quark operator of spin n is given by 



O 



Ml,M2— Mn . 



X 



,n-l S 



^(x)^ 1 D^ 2 ■ • ■ D^—^x) + trace terms 



In the singlet case there are two operators. The quark operator is represented by 



O 



1 

-I 

2 



n-l S 



■^(x)7 Ml D M2 • • • D lln ip(x) + trace terms 



and the gluon operator is given by 

1 



F^(x)D^ 2 ■ ■ ■ D^F^ix) + trace terms 



(2.1) 



(2.2) 



(2.3) 



In these composite operators ip and Fg v stand for the quark field and the gluon field tensor 
respectively. The in Eq. fl2.1|) represent the generators of the flavour group and the index 
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a in Eq. ( |2.3|) stands for the colour. Further the above operators are irreducible tensors 
with respect to the Lorentz group so that they have to be symmetric and traceless in all 
their Lorentz indices //j. From the operators above one can derive the Feynman rules for the 
operator vertices in the standard way (see e.g. ||, p~SH , |16|). This derivation is facilitated 
if the operators are multiplied by the source 

J = ^Ml^M2 ' ' ' ) (2-4) 

with A 2 = in order to eliminate the trace terms in Eqs. ( |2.1| )-( |IO| ). Hence all operator 
vertices in momentum space are multiplied by a factor (A ■ p) n . For the computation of the 
OME's denoted by 

Aij = 0'(P) I Oi I J(P)) (2-5) 

with i,j = q, g we choose the Feynman gauge except for the one-loop graphs for which we 
take the general covariant gauge. For this choice the gluon propagator equals 

A£(*O = ^(-*r + (l-0^). (2.6) 

The matrix element ( |2.5| ) has to be considered as a connected Green function with the 
external legs amputated but with the external self energies of the partons j included. In 
this paper all quarks and gluons are taken to be massless and the external momentum p is 
off-shell (p 2 < 0) in order to get finite expressions for the OME's. This choice implies that 
the OME's are not gauge invariant so that they cease to be ordinary S-matrix elements. 
Moreover they acquire unphysical parts which can be split into two classes. The first class 
originates from the fact that the equations of motion (EOM) do not apply anymore, which is 
the case for both non-singlet and singlet OME's. The second class can be traced back to the 
mixing between so-called gauge invariant (GI) or physical (PHYS) and non-gauge invariant 



(NGI) operators which originate from the Yang-Mills (here gluonic) sector (see 17], [fig] , [19 
p0|). Therefore this second class of unphysical operators only shows up in the singlet case. 
The first class does not contribute to the operator renormalization constants in contrast to 
the second class which affects the aforementioned constants via the mixing between the GI 
and NGI operators as we will see below. The way to deal with this mixing is described in 
T§, 0, and |Z2|. 



The calculation of the Feynman graphs corresponding to the physical operators of Eqs. 
(j2.1|)-(P75D, which are depicted in the figures in 0, proceeds in the standard way. The 



figures and definitions for the corresponding unphysical operators are given in fL5| . The 
Feynman integrals reveal ultraviolet divergences which are regularized using the method of 
N- dimensional regularization. In this way the above divergences show up in the form of 
pole terms of the type (l/e) k with e = N — 4. In ||- |J it was sufficient to evaluate the 
one-loop graphs up to finite terms and the two-loop graphs up to single pole terms in order 
to get the second order anomalous dimensions. Here we have to include terms proportional 
to e in the one-loop expressions and the two-loop graphs have to be computed up to terms 
which are finite in the limit s — > 0. The way to compute the two-loop Feynman integrals up 
to finite terms is presented in and in Appendix B of [23| . We used the program FORM 



24 to do the necessary algebra. 
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To check our results for the Feynman diagrams it is useful to have explicit expressions for 
the pole terms in e. Therefore we will now present the OME's expressed in renormalization 
group coefficients, which are defined in [[ll]] and |T2]. The explicit formulae for the physical 
and unphysical OME's can be found in the Appendices A and B respectively. Further it 
is implicitly understood that all quantities in the main text, in particular the anomalous 
dimensions 7^ = q,g), are Mellin transforms but to avoid additional indices we do not 
write a superscript n to indicate this explicitly. (Another way to interpret the formulae 
is that the OME's are given in parton momentum fraction (z) space when the anomalous 
dimensions are replaced by minus the corresponding Altarelli-Parisi splitting functions and 
the multiplications in moment space are replaced by convolutions.) We have written the 
OME's in such a way that all renormalization group coefficients appearing in the expressions 
below are renormalized in the MS-scheme. Up to order a 2 s the non-singlet and the singlet 
OME's can be decomposed as 



A 



qq 



r,PHYS 

qq 



A • p 

p2 



A 



r,EOM 

qq 



(A ■ p)> 



(2.7) 



where A r q '^ HYS and Ag' g EOM with r = NS, S stand for the physical and unphysical parts 
respectively. The latter enter due to the breakdown of the equations of motion. The non- 
singlet OME can now be expressed into renormalization group coefficients as follows 



.4 



NS,PHYS 

qq 



1 + a s S £ 



-P 



p W 2 T:!~ns,(o) 



NS,(1) _i_ r _NS,e,(l) 
u qq ' CLL qq 



) £ ^{i(7^ (0) ) 2 -^ (0) } 



el2 hq 



NS,(0) NS,(1) 

qq qq 



~d a NS ' (1) , 



NS,(2) 
l qq 



2/3 af- (1) 



"qq 



+7, 



NS,(0) a NS,e,(l) 



qq 



qq 



A NS, £ ,(1) . 
; u u qq 

^ — Tt — z t 



(2.8) 



where S e denotes the spherical factor given by 

-e 



exp 



7£ — ln47r 



(2.9) 



which originates from iV-dimensional regularization. In the expressions above all quantities 
that are unrenormalized with respect to operator, coupling constant a s and gauge constant 
£ renormalization are indicated by a hat. The finite terms are written in such a way that 
after all renormalizations the non logarithmic terms (with respect to ln(— p 2 /fJ?)) become 
equal to a^ s '^ (k = 1,2). Further we have introduced a shorthand notation for the strong 
coupling constant so that 
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An 



El 

Air 



(2.10) 



The coefficients (3q and which originate from coupling constant and gauge constant renor- 
malization, are given by 



a, = a. 



'l + a s S £ (2(3 -)~ 



(2.11) 



1 + a 



(2.12) 



where 



A) — 







~ 3 n f T f . 



(2.13) 



In QCD (SU(N)) the colour factors are given by C F = (N 2 - 1)/2N, C A = N, T f = 1/2 



(k) 

and rif stands for the number of light flavours. Finally the 7^ denote the coefficients of the 
order a k s +1 terms appearing in the series expansions of the anomalous dimensions. Using the 
same notation we can also express the unphysical part of the non-singlet OME in Q2.5p in 
the aforementioned renormalization group coefficients 



A 



NS,EOM 
QQ 



OisS, 



~P W2 



£ \ 2 
fi 2 



6 NS,(1) + £& NS,e,(l) 



qq 



qq 



+a 2 s S 2 ( 



11- 

^ &l s ' {1) 



)4i/ 7 NS.(o) 6 NS,(i) _ 2/3 Ml) 



QQ 



<- 



d &3 £ '« 



-4 + C (2) - 2 « s - (1) + 7 9 T (0) C £ ' (1) 



qq 



€=1 



(2.14) 



Since the singlet OME is only computed up to second order the unphysical part arises 
from the fact that the equations of motion are not satisfied. In this order there is no need 
to introduce NGI operators. The singlet OME can be split into nonsinglet (NS) and purely 
singlet (PS) parts 



A s 
qq 



A m + A 
qq qq 



PS 



(2.15) 



where the purely singlet physical OME is 



.4 



PS,PHYS 

qq 



arSlf— =-) 



I/I T (o) T (o)\ 

P 2 \ o Iqg Igq J 



if 1 



PS,(1) 



Iqq 



+ ^(0) a (l)\ +a PS,(2) , w , 
' Iqg a gq J ' u qq Iqq a 



(0)„£,(1) 

qg gq 



(2.16) 



and the purely singlet unphysical OME is 

-' r ^[I/ 7 (0) 6 (l)l + 6 PS,(2) + (0) 6 e,(l 



— i 2 

^PS,EOM _ ^C^t r \s 



The mixing with NGI operators does occur in A qg ^ v which can be written as 

A iPHYS T (l) , 4EOM T (2) , 4NGI T (3) 

^■qg,fiv si-qg ^ ^ v T ^- qg A fj,u ^ ^qg A fiu • 



In this expression the tensors are given by 



T (i) 



l + (-l) n 



5V 



A-p 



+ 



(A- P y 



(A- P y 



(2.17) 



(2.18) 



(2-19) 



and 



T (2) 
- 1 fiu 



PpA v + A^p v A^p 1 



P 



A-p 



(A- P y 



(a-pY 



T (3) 



A-p 



+ 2- 



(a-pY 



(A ■ pY 



For later purposes we also define the tensor 



T (4) 



1 + (-1)" [Pfl A u + 2 A^A^ 2 



A-p 



(A-pY 



(A ■ p)' 



The above tensors satisfy the following relations 



V 



T« = (< = 1,2), ^T«^0 (i = 3,4), 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



V 



yjJ = o (* = i,2,3), pYt£Vo 



(2.24) 



Using these relations one can show that the OME's in Eq. fl2.18Q satisfy the following 
Ward-identities (WI) for i = q 



p^A 



ig,/j,u 



~Pv + 



A u p 2 
A-p 



n A NGI 



(A • pTA\ 



KJ 



(2.25) 



p»p u A ig ^ = . 



(2.26) 



The WI in Eq. ( j2.25| ) shows that the unphysical part in Eq. ( 2.18|) , given by A q , , is due to 
the NGI operator O a in Eq. ( |B1|) . The second term in Eq. ( |2.18| ), given by A^ M , is also 
unphysical due to the fact that we cannot apply the equations of motion. The physical part 
of the OME in Eq. ( gTp becomes 
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/jPHYS 
A <?9 



l qg 



2 c ,2, P 2 \e\ 1 f 1 f„,(0) o/ (0) 



, [MV 7 (o) 7 (o) + (o) 7 (oA_/3 o7 (o)j 
• K ^2 ' l e 2\2\ 1( N J i9 ^ ! qg '99 J l iJ ! qg J 



e^2 q9 qg qg 99 w * eg 5 J 

+a£> - 2/^f + 7 £>a# 1 > + 7 £>a£< 1 > 



The unphysical parts are given by 



(1* 



+ aZS 2 (—) £ \-W% {l) + 7 (0) 6 (1) - 25 n 6 (1) 
/i 2 99 ! qq qg ^^v u qg 



and 



/i 

W 0) aM 1} 



(2.27) 



^2.28) 



(2.29) 



The coefficients with a subscript A originate from the NGI operator O a whose matrix element 
will be presented below when we discuss the renormalization. 
The next OME A gq can be decomposed as 



A gq ={4A p ™ + t^PA*™](A. P r, 



(2.30) 



where the last term again represents the unphysical part due to the breakdown of the 
equations of motion. The physical part equals 



jPHYS 

gq 



= a s S £K 

fi 2 



~(0) i fl (l) _J_ £a e,(l) 

'gq ^ a gq ^ tu gq 
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1 rl 



£ ^2 



+7»( 



(0) / 6,(1) 



,(2) 



2/^f + ^^ + ^f 



+ a 



Bg / 



•H 



(2.31) 



The coefficients with the subscripts A and B refer to the NGI operators Oa and Ob presented 
in Eq. ( pi|) and Eq. ( pB5|) respectively. The latter shows up for the first time in order a 2 
in the physical part of A gq . The corresponding OME's will be given below when we discuss 
the renormalization of the physical operators. The unphysical part of Eq. ( [2.30 ) equals 



A 



EOM 
99 



A(-fr /2 [f«+^f 



+aiS, 



2c2/ P 



di > 99 99 99 99 99 99 

~d 6 £ ' (] 



(2.32) 



The last OME A ggtfiv has a similar decomposition to the one presented in Eq. (|2.18| ) 

A 4PHYS T (1) , 4EOM T (2) , 4NGI T (3) 



(2.33) 

In particular it satisfies the same Ward-identities as listed in Eq. (|2.25|) and Eq. (|2.26|) for 
% = g. Hence the second (unphysical) part in this expression originates from the breakdown 
of the equations of motion whereas the third term again originates from the NGI operator 
Oa in Eq. (Bl). The physical part of the OME is equal to 



.4 



PHYS 
99 



i+^(^r /2 [^+aS+^ (1) 



+a 2 S 2 ( 



I(if( 7 (0)^2 



'99 



+ ~(0)^(0) 
'91 '99 



£ ^2 
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.</.</ 



The unphysical parts are given by 



d£ 



^eom _ Q s g £ ^ ^ 5(1) _|_ gfo^ 1 



99 



and 



i{->s?ts/-wa?+->a(»2+e 



• d &<!> 



(2.34) 



^2.35) 



n 99 ~ as ^^^2 



~P ,e/2 



LJ® 4. 



9^ 



^2 



^<i a 



(i) 



£ 12 9 ^ _ l ^ a gA + %9 a gA + 7<,Al a Ayl + a uMJ ~ <5 ^ ^/ 



.(2) 
l 9 A 



^P^gA +Tgg a gA + lgA\ a AA + a uA , 



d a £ ' {1) 
a a gA 



^2.36) 



respectively. In these expressions we observe that there are quantities with the subscript A 
and with the subscript to. The latter originate from the NGI ghost operator presented 
in Eq. flB2"f ). The renormalization of the above OME's involves the mixing of the NGI 
(non gauge invariant) operators Oa, Ob and O u mentioned above with the physical (gauge 
invariant) operators in Eqs. (|2.1| )- (|2.3| ). Therefore we have to compute the matrix elements 
in Eq. (|2.5|) where the physical operators indicated by i = q, g are replaced by the NGI 
ones labelled by i = A, B, u. To get the physical OME's which are finite up to order a 2 s the 
unphysical ones have fortunately only to be calculated up to order a s . 

The operator renormalization proceeds as follows. First we have to perform coupling con- 
stant and gauge constant renormalization. This is achieved by replacing the bare constants 
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by the renormalized ones by substituting Eqs. (|2.11|) and ( |2.12| ) in the above expressions for 
the OME's. Subsequently the OME's have to be multiplied by the operator renormalization 
constants to remove the remaining ultraviolet divergences. The most simple case is the 
renormalization of the non-singlet OME's since here we do not have mixing with physical 
and NGI operators. This renormalization is achieved by 



4 NS,i 
^qq 



(z-x s K s ^ ( 2 - 37 ) 

with i = PHYS and EOM. The inverse of the operator renormalization constant equals 

1 



-1\NS 
Iqq 



_^ N S,(0) 



NS,(0h2 _ B m,(o)\ _ _L NS,(1) 

£n 2 uw qq J 2e q 



(2.38) 



In the case of the singlet operators A^ q l there is no mixing between GI and NGI operators 
(at least up to order af). The singlet OME's are finite for i = PHYS and EOM when 



.4 



qq 



qq- ^qq 



-1\ Ai 
qg-^gq > 



^2.39) 



and the renormalization involves the matrix Zij due to mixing with unphysical operators. 
Hence we have to invert this matrix in order to perform the operator renormalization. The 
first renormalization constant in Eq. (|2.39|) can be split as follows 



-1\S 



1\NS _|_ /^-InPS 
qq \~ Iqq \ I qq 



where the first part is equal to the expression ( j2.38j ) and the last part equals 



(Z 



-1\PS 
Iqq 



1 rl 



\o Iqg Igq J 



PS,(1) 



2s 



Iqq 



(2.40) 



(2.41) 



The constant (Z l ) qg in Eq. ( 2.39Q also shows up in the renormalization of the next OME, 
namely 



A' 



(Z 



-1\S 



Ai n + (Z 



-1\ 



A 1 



^2.42) 



for i = PHYS, EOM and NGI. It can be decomposed into renormalization group coefficients 
as follows 



(Z 



><l!l 



a,S f 



- ±~(o) 
e Jg9 



rl _/I 7 (o)( 7 (o) + 7 (oh_« 7 (o)\ 
2 \ <-> iqg V 'qq ~ hg > 1 " <qg / 



__ 7 W 
2e Jq9 



For the renormalization of A™ YS in Eq. ( |2.31| ) and v4^ M in Eq. 
contributions of the NGI OME's. Here the finite expressions are given by 



4q = (Z^gg^q + v{A\ q + & Bq )} + [Z^) gq A 



qq 



(2.43) 

we also need the 
(2.44) 
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where i] is defined as a source multiplying the NGI operators in the effective action (see Eq. 



(2.8) in ||15|| ). Here we will write this quantity as follows 

i] = a s S t 



1 (0) 



2 C2 T 1 / f (0) , (0) (0) 



an 



(2.45) 



If we sandwich the operator Oa between quark states (see Eq. (f2.5|)) we get 



(o) , (i) , 

lAn + & An + ? a 



Aq 



and 



5 (i) 



eb 



,£,(1) 



^2.46) 



(2.47) 



up to lowest order. We proceed in a similar way for the operator Ob and the OME's read 

(2.48) 



1 (0) , (1) , e,(l) 



and 



/(HUM — p. Q 
A Bq - a sOe\ 2 



P\e/2 



eb 



£,(1) 



^2.49) 



Notice that A Aq and can be decomposed in a similar way to A gq in Eq. (|2.30|) . The 
operator renormalization constant appearing in the last part of Eq. ( 2.44|) can now be 
written as 



(Z 



-In 



I gn 



a,S f 



~2e lgq 



- ±~(o) 



(i) 



+ «^^i{k ) (7S ) + 7S ) )-/5o7S ) } 



-e 2 ^2 



^2.50) 



The first constant [Z 1 ) g9 in Eq. ( |2.44| ) also shows up in the renormalization of the OME 



A gg = {z- x ) m A\ g + (z-% [A gg + n (A\ g + A; 



(2.51) 



for i = PHYS, EOM and NGI. If we express this constant into the renormalization group 
coefficients in the following way 



[z- x ) 



nn 



1 + a s S £ 



1 



199 



1(I( 7 (0))2 + I (0) (0) _ g (0)1 _ 1(1) 



(2.52) 
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we obtain the finite physical (PHYS) and unphysical parts (EOM, NGI) of expression ( [2.511) . 
In order to obtain these finite parts one has first to calculate the OME's which emerge when 
the NGI operators Oa and are sandwiched between gluon states (see Eq. ( |2.5| )). They 
can be decomposed as follows 



A _ 4 PHYS T (l) , 4EOM T (2) , 2NGI T (3) , JWI T (4) 



(2.53) 



with 



< YS = «A(-j-r/ 2 



,7a 9 + a A 9 + £a Ag 



(2.54) 



A Ag — a * D < 



~" A /I 2 ] 



1.(1) , Fh e,(l) 
°Ag + £0 Ag 



[2.55) 



A Ag — 1 + a s S l 



-V \e/2 



fi 2 



- (0) (1) £,(1) 

+ a AA + £a AA 



[2.56) 



and 



A 9 - "s^eV ^ 2 



r (l) , p ^,(1) 
C A S + £C Ag 



[2.57) 



Similar algebraic expressions are obtained for O u . They are given by 



A 4PHYS T (1) , 4EOM T (2) , 2NGI T (3) , % WI T (4) 

■"-uig^u -"-ujg ± fiv ~ ^uog ± \iv ~ ug fxv ~ ^ug ± fiu 



[2.58) 



with 



iPHYS _ ~ o 

A wg ~ a s^ey 2 

A* 



up 



^2.59) 



iEOM _ f. c 
■^ujg u>s^e\ 2 



+ e6E ,(l) 

u ujg ' ug 



[2.m) 



(2.61) 



and 



4 wi _ - q 



Ae/2 



^2.62) 



Notice that neither Aa 9 ^ v in Eq. (|2.54j ) nor A ug ^ v in Eq. ( |2.59| ) satisfy the Ward-identities 
(WI) in Eq. ( |2.25| ) and Eq. ( |2.26| ). However if we add them according to Eq. fl2.51| ) the 
terms proportional to the tensor Tj£) cancel and the Ward-identies are restored. 
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From the residues of the single pole terms appearing in the physical operator renor- 
malization constants Zij with i,j = q,g given above one can now read off the anomalous 
dimensions presented in the MS-scheme. Here we agree with the results published in || (see 
also ||15|| ). Therefore we have a check that the residues of the single and double pole terms 
are correct. By comparison with the algebraic expressions given above one can obtain all 
renormalization group coefficients like a\j , etc.. Explicit expressions for the physical 
and unphysical OME's can be found in Appendices A and B respectively. They contain all 
the finite second order terms which survive in the limit e — > 0. These terms have not been 
calculated previously. 

Note that a remarkable property is found for the non-singlet OME A^ q s ' PHYS in Eq. ( |2.8| ) 
for which the full expression is given in Eq. ( |A3|) . The first moment of this OME is equal to 
unity up to second order in a s provided we choose the Feynman gauge (£ = 1). This result 
is expected for the on-shell expression as it is a check of the Adler sum rule |25] . However 
the fact that all the coefficients of the terms in (— p 2 //i 2 ) are zero shows that the sum rule is 
true up to order a 2 s for the off-shell expression which is not an S-matrix element. Finally we 
want to comment on the use of the unrenormalized expressions given in Appendix A for the 
computation of the three-loop OME's. Since the external quark and gluon legs are off-shell 
one can insert our results into the Feynman integrals for one-loop graphs. In this way one 
gets expressions corresponding to some of the three-loop graphs. It is however clear that 
the most difficult Feynman integrals belonging to the non-planar diagrams, where all quark 
or gluon lines cross over, remain to be done. 
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APPENDIX A: 



In this Appendix we present full expressions for the two-loop corrected operator matrix 
elements computed from the Feynman diagrams depicted in ||. The second order con- 
tributions are calculated up to finite terms which survive in the limit e — > 0. The OME's 
presented here are unrenormalized and external self-energy corrections are included. In these 
expressions definitions of the Riemann zeta-f unctions ((n) and the polylogarithms Li n (z), 
Sn,m( z ) can be found in P5| . Also the distributions (1/(1 — z)) + and (ln(l — z)/{l — z)) + 
are written simply as 1/(1 — 2;) and ln(l — z)/(l — z) respectively to shorten the formulae. 
Note that the OME's given in the text are the moments of the functions listed here so 

A§ = £ dzz^A^z,^^), (Al) 

where for simplicity we have not written the moment index n on the functions. Also to 
simplify the expressions we define the phase-space factor 

F = %S A ^. (A2) 
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We first split into physical and unphysical parts following the notation in Eq. 
The physical part is 



.4 



NS.PHYS / 
<?<? 



(---£;)=*< 1 -*> 



/i 2 e 



+F C F [-{ - 4 - 4z + — — + 65(1 - z)\ + 4 - 6z 



+2(: 



1 - z) [In z + ln(l - z)\ - 5(1 - z) [7 - 4£(2)] + 

2 



1-2 

*)[C(2) 

+{lnz + ln(l - z)} 2 ] + 5(1 - z)[7 - ~C(2) - 4C(3)] + 1(1 - |) 



+ £ {l - 2z + (2 - 3z)[ln z + ln(l - z)] + 1(- 



x[2 



+F 2 



1 ln(l - z) 



lnz 



1 - z 



1 - z 



+ — + 5(l-z){-2 + C(2)}]} 



-ACI [40 + 8z - - 25(1 - z)[9 - 16C(2)] 
1 1 — z 



81nz(3 + 3z 



+321n(l-z)(l + z-- ; 

44 2 
+CUC J ,[— (1 + z - ) - 225(1 - z)] 

+n f C F T f [-f(l + z- T ^-] 

+ -{C 2 [-40 - 4z + + 5(1 - z)[^ - 36C(2) - 8C(3)] 

£ L 1 — Z 2 



55(1 -z)]} 



-41n(l - z)(l - llz + 



1 -z 



+4 lnz(9 - 5z - 
-21n 2 z(7 + 7z 



9 



+ 241n 2 (l-z)(l + z- 



1 



l-z' 



16 lnz ln(l — z) 

1 — z 

-8Li 2 (l-z)(l + z)] 



1 - z y 

^ ^ r 158 22 238 1 r/1 , r 325 44 A . . 
-CUC F [ — — + —z + — + 5(1 -z)[—- — C(2) 



9 9 
■12C(3)]+4(l + z 



9 l-z v /L 6 
T ^)[C(2) + ^ln(l-z)-lln 2 z] 



, ,34 34 44 1 

+ lnz( 1 z ) 

v 3 3 3 1 - z n 

_ ^ r 88 56 56 1 r/1 u 58 16., nVl 
+ n/ C F T f [- - -z - -— + 6(1 - ,)[_ + T C(2)] 



--{lnz + 21n(l-z)}(l + z 



1 - z 



)]} 



188 4 

h -; 

3 3 



+C|[^ + ^z 



56 , w 32 „ 4 v 

_ + C (2)(20-4, + T ^- — 



1 541 Q7 

+24C(3)(1 - z + 2z 2 - T --) + 5(1 - z)[- a — + ^C(2) 
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74 98 
+54C(3) - ^-C(2) 2 ] + ln(l - 2) (-38 + 382 + — 

O 1 z 

+81n(l - ^)C(2)(1 + 2 - 22 2 - — *— ) - ln 2 (l - z)(13 - 31z 

1 — z 

+-^) + ^ln 3 (l-2)(l + 2 2 



l-2 y 3 v ' y \-z' 
-41n(l - z)Li 2 (l - z)(3 - 2 + 4z 2 ) - lnz(^ + ^z 11 



3 3 \-z' 
10 

+ ln<(2)(2 - 14.2 + I62 2 ) - 21nzln(l - z)(3 - 17 z + ) 

1 — z 

+2hi2ln 2 (l -z){3 + 3z- -^-) + In 2 2(11 - 232 - ^-z 2 

1 — z 3 

l< ' ln 2 2ln(l- 2)(6 -22 + 82 2 + — ^)-ln 3 2(5 + 52 



\-z' v l-2 y 

X 6 1 2 

-) -8 In z\A 2 {l-z) (1 - 32 + 42 2 + ) 



3\-z' " ' v /v l-2 y 
-4Li 2 (l - z)(6 - 122 + y~) + 4Li 3( X - - 7z + 122 2 ) 
8 1 

+-(9 + - + 12z + 42 2 )[hi2ln(l + z) + Li 2 (-*)] - 16[lnzLi 2 (-z) 

3 2 

-2Li 3 (-z)](l - ) - 8S 12 (1 -z)(l-7z + Qz 2 + -i-)] 

i. Z L Z 

1 670 1 98 

+C A C F [-[941 - 5802 - - -C(2)(l - 17* + I62 2 + 

-2C(3)(10 + Az + 12, 2 - ^ + 5(1 - *)[-^ + ^C(2) 

4Q 2 11Q 

+28C(3) + ^C(2) 2 ] - 5 ln(l - z)(61 + 16z - -^-) + 2 ln(l - z)C(2) 
5 9 1 — 2 

5 22 2 

x(l + 32 + 42 2 - ) + {— ln 2 (l -z)- In 3 z}(l + z- - ) 

1 — z 3 1 — z 

5 1 22 
-2 In z((2) (3 + z + Az 2 - - ) + - In 2 2(47 + 47 z + 16z 2 - - ) 

1 — 20 1 — z 

2 19 1 101 
+-lnzln(l -z)(U + 5z- - ) + - \nz(U + 1762 + ) 

3 1 — z 9 1 — z 

+4hi2Li 2 (l -z)(2 + Az 2 - -^—) - 4Li 2 (l - z)(3z - — *— ) 

1 — z 1 — z 

-4(3 + + Az + \z 2 ) [In 2 ln(l + 2) + U 2 {-z)\ 

OZ O 

+ (1 - z + Az 2 - — !— ) [In 2 z ln(l - z) + 2 ln(l - z)Li 2 (l - 2) 
1 — 2 

-6Li 3 (l - 2)] + 8[ln2Li 2 (-2) - 2Li 3 (-2)](l - — *— ) 

1 — 2 

+S 12 (1 - 2)(14 - 62 + 242 2 - -i^)] 
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4 32 569 46 

+n / C i? T / [-[-22 + 5z + — ] + 5(1 - z)[— - -C(2) 

-8C(3)] + l{21n(l - z) + In 4(11 - 7z - 

[2C(2) + {2 ln(l - z) + In z} 2 ] (1 + z - ^)] 



-2(-l) n (C 2 F - \c A C F ){^{ - 8 + 8z - 41nz(l + z) 

+2(-l + ^ + —— ) [2C(2) + 4 In z ln(l + 2) - In 2 z + 4Li 2 (-z)] } 

1 — h 2- 



41 



-(1 -z)-8 ln(l - z)(l - z) + 2C(2)(3 - 3z + -z 2 ) 



+8(z + 



l + z 



) ln(l + z)[({2) + lnzm(l + 2) + 2Li 2 (-z)] 



-4(1 + z)[lnzln(l - 2) + Li 2 (l - z) + Li 3 (-z)] 

1 8 
— In 2(37 + 25^) - In 2 z(6 + -z 2 ) - 2 In 2 z ln(l + z) (3 - z - 



l + z' 



-2C(3)(3 + z- -1— ) + ^(- + 3z + 4z 2 ) [In z ln(l + z) + Li 2 (-z)] 
1 + z 3 z 



-8S 



121 



-z)(-l + 3z + 



l + z- 



+ 4(l-z 



1 + z /L 4 



)h^z 



-lnzln 2 (l - z) -21nzln(l - z) ln(l + z) + -ln 2 zln(l - z) 

-2 ln(l - z)Li 2 (l - z) - 2 ln(l + z)Li 2 (l - z) + lnzLi 2 (l - z) 

- In zLi 2 (-z) + 214(1 - z) + S 12 (l - z) + S 12 (z 2 )])] . (A3) 

Here the factor (—1)™ originates from the non-planar diagrams (namely k and o in figure 2 
in the singlet paper of 0). It multiplies that part of the matrix element which is needed for 
the mass factorization of physical processes with two identical quarks in the final state. 
The unphysical part is given by (see Eq. (|2.14| )) 



^NS.EOM 



qq 



fi' 2 : 



F C F 



4z - 2(1 - i) + e[2z - (1 - £)}{l + lnz + ln(l - z)} 



1 

+F 2 [-{C 2 F [-lQ-8z+ 16z(-21n(l - z) +lnz)] 



+F 2 










. .20 


+C A C F [— - 




32 


+c 2 F [ 





88 16 32 

yz] + n/C>2>[-y + —z 



-« + C(2)(- 



32 , 
16z z 2 



44zln(l - z) 



-24z ln 2 (l - z) - - In z(10 - 17 z) + 8(1 - 4z) In z ln(l - z) 
3 

16 1 16 

+— (- -3z-2z 2 )[lnzln(l + z) + Li 2 (-z)] +ln 2 z(12z + — z 2 ) 

O Z O 

+16Li 2 (l - z)(l - 3z) + 8z(l - z){-2 ln(l - z)C(2) + 2 In z((2) 
- In 2 z ln(l - z) - 4 In zLi 2 (l - z) - 2 ln(l - z)Li 2 (l - z) 



17 



+6C(3) + 6Li 3 (l - z) - 6S 12 (1 - z)}] 

26 232 i 
+C A C F [-— + —z + 4C(2)(1 - 2z + -z 2 ) + lnz(6 - 28z) 



16 



ln(l - z)(l + 4z) - 4[lnzln(l - z) + 2Li 2 (l - z)](l - 2z) 



1 



--z 2 In 2 z- ^ -3z-2z 2 )[lnzln(l + z) + Li 2 (-z)] 

+82(1 - z){- In 2 zln(l - z) + C(2)[ln(l - z) - lnz] + 2 lnzLi 2 (l - z) 
-3C(3) + ln(l - z)Li 2 (l - z) - 3Li 3 (l - z) + 3S 12 (1 - z)}] 
+n / C F T / [-^(l + 4z) - |(1 - 2z){lnz + 21n(l - z)}] 

-y (-l) n (^l - ^AC F ){l - z + ^ 2 [-2C(2) + In 2 z] 

- lnz(l - 2z) + (- - 3z - 2z 2 )[lnzln(l + z) + Li 2 (-z)]} 

The purely singlet (PS) OME is split in the same way as the non-singlet one in (]2] 
the physical part we have (see Eq. (j2.16| )) 



A 



PS,PHYS 
11 



( h 



/z 2 : 



rlfi r 4 4 1 

F 2 n f C F T f [— [1 + — - z - -z 2 + 2lnz(l + z)} 



+ -f-(16 + --13z-8z 2 ) + 21n(l-z)(l ' 



32 



4 ^ 

z z 

3 ; 



+ In z(13 + — + 13z) + 4(1 + z) [In z ln(l - z) + - In 2 z + Li 2 (l - z)] } 

316 20 256 , r , . , 4 4 9x 

+ — + - - — z - 40z 2 + 4{C(2) + 2 ln(l - z)}(l + — - z - -z 2 ) 

16 , . , . 5 9 , 9 . 16 8 9x 

+ — ln(l - z)(16 + - - 13z - 8z 2 ) + 2 In 2 z(37 + — + 41z + -z 2 ) 

8 8 10 

+81nzln(l - z)(13 + — + 13z) + -lnz(67 + — + 50z) 

+ 16Li 2 (l - z)(6 + — + 7z + \z 2 ) + 16(1 + z)[- lnz((2) 

3 

+2 ln(l - z)Li 2 (l - z) + In z ln 2 (l - z + - In 2 z ln(l - z) 

+ L in 3 z + 2 In zLi 2 (l - z) - 2Li 3 (l - z) + Si 2 (l - z)] 

and the unphysical part (see Eq. ( |2.17| )) becomes equal to 

-p 2 1- 



.4 



PS,EOM 



[l 2 E' 

-32 



F 2 n/C^T/J— { - 2 + z + z 2 -lnz(l + 2z)} 
-72 + 32z + 40z 2 + 32 ln(l - z)(-2 + z + z 2 ) - 16 In z(6 + 7z) 



(A4) 
. For 



(A5) 
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-32(1 + 2z)[lnzrn(l - z) + - In 2 z + Li 2 (l - z)\ 



(A6) 



The next OME is split into three pieces according to Eq. ( |2.18|) . The physical part is given 
by (see Eq. ( gg7p ) 



4PHYS/ 
A ig ^' ,,2 



p 2 1 



F n f T f [-{l - 2z + 2z 2 } + 4 + 4(1 - 2z + 2z 2 )[lnz + ln(l - 2)] 
+£{2[lnz + ln(l - z)] + {l-2z + 2z 2 )[C(2) + {lnz + ln(l - z)} 
+F 2 n f [\ [C A T f [^(U + - + 2z - 9z 2 ) + 32 lnz(l + Az) 

& o Z 

64 

+32 ln(l - z)(l - 2z + 2z 2 )} - n f T 2 — [l - 2z + 2z 2 } 

3 

-8C F T f [1 - Az - 4 ln(l - z) (1 - 2z + 2z 2 ) + 2 In z(l - 2z + 4z 2 )] } 
+ -{C A ^>[|(26 + — - 79^ - 20z 2 ) + 241n 2 (l - z)(l - 2z + 2z 2 ) 



+32C(2)(1 - z + 2z 2 ) + - ln(l - z)(17 + - + 38z - 52z 2 ) + 32(1 + Az) 

o Z 

x In z ln(l - z) + 8 In 2 z(3 + lOz) + — In z(25 + - + 13z + 14z 2 ) 

3 z 

+64z(3 - z)Li 2 (l - z) + 16(1 + 2z + 2z 2 )[lnzln(l + z) + Li 2 (-z)]] 

+C F T f [-32 + 36z + 8(1 - 2z + 2z 2 )[-2((2) + 3 ln 2 (l - z)] 

-4 In z(5 + 4z + 8z 2 ) + 8 ln(l - z)(3 + Az 2 ) - 12 In 2 z(l - 2z + 4z 2 ) 

+32 lnz ln(l - z)(l - 2z + z 2 ) + 16Li 2 (l - z)(3 - 6z + 4z 2 )] 
32 1 

+n / T / 2 y[-(2 - lOz + 10z 2 ) - (1 - 2z + 2z 2 ){lnz + ln(l - z) 

„ m r 2900 20 136 2492 2 > . w 16 172 
+C A T f [^ r + — - —2 - + C(2)(12 + — + 16z 



27 



27 



27 



3z 



28 



+16z 3 ) + 8C(3)(1 - z - 4z 2 + 6z 3 ) + — ln 3 (l - z)(l - 2z + 2z 2 

3 

2 16 
+ - ln 2 (l - z)(23 + — + HOz - 138z 2 ) + 81n(l - z)Li 2 (l - z)(4 
3 



z 



+23z + 2z 2 - 2z 3 ) + - ln(l - z)(68 



60 



283z + 94z 2 



+8 ln(l - z)((2)(3 - 7z + 8z 2 - 2z 3 ) + 8 lnzln 2 (l - z)(3 + lOz + 2z 2 ) 

+ - In z ln(l - z) (83 + — + 86z + 22z 2 ) + 4 In 2 z ln(l - z) (5 + 21z 

3 z 

-2z 3 ) + - In 3 z(7 + 22z) + - In 2 z(61 + - + 49z + 51z 2 - 6z 3 ) 
3 3 z 

8 30 

+16 In zLi 2 (l - z)(l + 7z - 2z 3 ) + - In z(157 + — + 127z + 59z 2 ) 

9 z 



+81nzC(2)(2 + 3z + 2z 3 ) + 16Li 2 (l - z)(5 



3z 



+ 4z + 



37 



19 



-8(1 + 4z + z 2 - 2z 3 ) [In z ln(l + z) + Li 2 (-z)] 
-8Li 3 (l - z)(8 + 13z + Uz 2 - 6z 3 ) + 8S 12 (1 - z)(5 - 5z + 12z 2 
-6z 3 ) + 8(1 + 2^ + 2z 2 ){2 ln(l + z)Li 2 (l - 2) + 41n(l + z)U 2 (-z) 
+21n(l + z)C(2) + 21nzln(l - z) ln(l + 2) + 21nzln 2 (l + 2) 



+- In 2 z ln(l + z) + In zLi 2 (-z) - Li 3 (-z) + 6S 



12 



Si 2 (^ 2 )}] 



+CV7)[-34 + 58z - 2C(2)(9 - 12z + 8z 2 ) 

+4(1 - 2z + 2z 2 )[-10C(3) + - ln 3 (l - z) - 2 ln(l - z)£{2)] 

+41n 2 (l - z)(5 - 2z + 6z 2 ) + 81n(l - z)Li 2 (l - z)(5 - lOz + 6z 2 ) 

-16z(l - 3z) ln(l - z) + 41nzln 2 (l - z)(7 - Uz + 10z 2 ) 

+4 In z ln(l - z) (7 - 16z + Az 2 ) + 4 In 2 z ln(l - z) (5 - lOz + 4z 2 ) 

-y In 3 z(l - 2z + 4z 2 ) - In 2 z(13 + 20z + 24z 2 ) 

+321nzLi 2 (l - z)(2 - 4z + 3z 2 ) - 21nz(34 - 15z + 24z 2 ) 
-41n<(2)(5 - lOz + 12z 2 ) + 8(3 - 6z + 2z 2 )[2Li 2 (l - z) 
-Li 3 (l - z)\ + 16Sia(l - z)(3 - 6z + 5z 2 )] 

-n/T?f [-(13 - 56z + 56z 2 ) - -(1 - 5z + 5z 2 )[lnz + ln(l - z)] 
3 9 3 

+ (1 - 2z + 2z){ln z + ln(l - z)} 2 



The unphysical part (EOM) (see Eq. ( |2.28|) ) becomes equal to 



.4 



EOM 
19 



( h 



F n f Tf - 16z(l 



8ez(l-z)[lnz + ln(l - z)] 



1 1 28 

+F 2 nf [-{c A T f [4& + —z 



272 



+ 161nz(l + 8z + 8z 2 )] 



+C F T / [-32 + 32z - 1282(1 - z) ln(l - z) + 64z(l - 2z) lnz] 
+n f Tj™z(l-z)} 

+C A T f [68 -™z- ™z 2 + 32C(2)(-z + z 2 + z 3 ) 



+16C(3)(z - llz 2 + 6z 3 ) + 16 ln(l - z)(3 + 



25 



34 



z 2 ) 



3 3 

+ 161n(l -z)C(2)(z + z 2 -2z 3 ) + 161nzln(l -z)(l + lOz + 6z 2 
+4 In 2 z(3 + 24z + 24z 2 - 4z 3 ) - 32z(l - z) lnzLi 2 (-z) 



50 



34 



+81nz(7 + —z + —z 2 ) + 16Li 2 (l - z)(l + 12 z + 4z 2 ) 
3 3 

+32(-l + 2z 2 + z 3 )[lnzln(l + z) + Li 2 (-z)] + 16(z - 3z 2 + 2z 3 ) 
x[- ln(l - z)Li 2 (l - z) - iln 2 zln(l - z) -21nzLi 2 (l - z) 



+ lnz((2) + 3Li 3 (l - z)] - 16Si 2 (l - z)(5z - llz 2 + 6z 3 



20 



+32z(l - 3z)Li 3 (-z) + 32z(l + z){-~ In 2 z ln(l + z) 

+ lnzln 2 (l + z) +ln(l + z)((2) + 2ln(l + z)U 2 (-z) + 2S 12 (-z)}] 
+C F T f [-16(l -z) + 32z(l - z)[2C(2) - 3 ln 2 (l - z)] 
-32 ln(l - z)(l + 3z - 4z 2 ) - 64z(2 - z) In z ln(l - z) 

+48z(l - 2z) In 2 z - 161nz(l - z + 8z 2 ) - 64z(3 - 2z)Li 2 (l - z)] 

64 5 
+n f Tf—z(l - z)[-- + lnz + ln(l - z)] 



(A8) 



The part which is due to the contribution of the non-gauge invariant (NGI) operator Oa 
(see Eq. ( p.29|) ) is equal to 



A 19 \ A 



-p 2 1 



9 ' 

H z e 



* - - 15/2. 

1 e 2 ^3z 



F 2 n f C A T 



— z 



-z 2 + lnz} 



+ -{ - |(7 + - - 13^ + ^ 2 ) - l 4 ln(l - z)(- - 3z + z 2 ) 
e L 3 z 3 z 

-161nzln(l -z)- 121n 2 z - 81nz(5 + — + z) - 16Li 2 (l - z 

3z 

104 10 122 „ 4,2 9 , r>/ , 

-— - - + —z + 4z 2 - -(- - 3z + z 2 )[C(2) + 2 ln 2 (l - z)] 

-16 ln(l - z)Li 2 (l -z)-~ ln(l - z)(7 + - - 13z + z 2 ) 

-81nzln 2 (l - z) -81nzln(l - z)(5 + — + z) - 121n 2 zln(l - z) 

3^ 



— In 3 z - 2 In 2 z(15 + — + 5z - -z 2 
3 v 3z 3 y 



161nzLi 2 (l -z) 



-| lnz(23 + - + lOz) - 4Li 2 (l - z)(10 + — + 4z - |z 2 ) 
3 z 3z 3 

-41n<(2) + 16Li 3 (l - z) - 8Si 2 (l - z) 



(A9) 



The next OME is split into physical and unphysical parts according to Eq. (|2.30|) . The 
former (see Eq. ( |2.31| )) is 



A 



PHYS 

9<1 



-p 2 1 

9 ' 

/r e 



F C F 



r4 r 2 1 4 4 
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The gluonic OME can be split into three parts according to Eq. ( |2.33| ). The physical part 
is equal to (see Eq. ( [2.34D ) 
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The unphysical part is given by (see Eq. ( |2.35| )) 
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The part which is due to the contribution of the NGI operators Oa and (see Eq. (|2.36|) ) 
is equal to 
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APPENDIX B: 



In this Appendix we present the non-gauge invariant (NGI) operators with their corre- 
sponding operator matrix elements (OME's). They are needed for the renormalization of 
the physical operators due to the mixing between them discussed in section 2. Two of them 
are given in Eqs. (2.9) and (2.10) of |15[] (for their construction see [0). These operators, 



which already show up in the case where only the gluonic operator in Eq. Q2.3|) is present, 
are given by 
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where £ and uo are the ghost and antighost respectively. In the expressions above 77 is defined 
in Eq. (|2.45 ). Further K{ and rji are calculated in |1| and are given by 



* = 5<-i)' + f- 



(n-2)! 



3 (n-2)! 
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and 
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respectively. If we also include the quark singlet operators in Eq. (|2.2j ) then we have to add 



U B 



x) = i n - Y S g^ k {x) 1 ^{T a ) kl A^{x)d^ ■ --d^tix) + 0{g i 
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Notice that the above operators are not BRS-exact in the strict sense of |T7| (see also [pTfl). 
This might affect the non-logarithmic terms in the renormalized OME's which we do not 
need here. Further the operators in the above equations are corrected up to order g 2 . This 
is sufficient to get finite two-loop OME's. However in order to carry out the renormalization 
on the three- loop level one has to compute higher order corrections to Eqs. ( |Bf| ), ( [B"2| ) and 
([BSD. The operator vertices corresponding to Oa in Eq. (|Bf| ) and Ob in Eq. (p2|) , albeit 
sandwiched between gluon states, are presented in Appendix A. 3 of |TJ|. In section 2 we 
need the Ob operator vertex when it is sandwiched between quark states. In this case we 
have to compute the quark-quark-gluon vertex which is given by 



= 0(T a )„4A"(A ■ k) 



n-2 



(B6) 



where k stands for the momentum of the gluon. 

We now list the OME's needed for the renormalization of the physical operators in section 
2 (for the notation see Appendix A). The two OME's referring to the NGI operators Oa 
and Os, when sandwiched between quarks states, can be split according to Eq. ( |2.7| ). The 
expressions in ( |2.46|) and ( |2.47|) become 
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respectively. Likewise we obtain for ( [2.480 and (|2. 19 ) 
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' 9 ' 

H z e 
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respectively. Next we sandwich the NGI operators Oa and O w between gluon states. De- 
composing the OME's according to (|2.53|) we get for Aa 9 the following results. The physical 
and unphysical parts (see Eq. ( |2.54j ) and ( |2.55| )) become equal to 
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5 1 
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respectively. The piece coming from the NGI operator Oa itself and the part due to the 
violation of the Ward-identity (WI) (see Eqs. ( |2.56|) , ( |2.57| )) are given by 
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respectively. 

The last OME due to can be also split into four pieces according to Eq. (|2.58|) . For Eq. 
( [2.59| ) we get 



-p 2 1 

/I 2 ' £ 



F C A z(l - z) - + ln(l -z)+\nz 



+£ i [{ln(l-z)+ln2} 2 + C(2)] 



The expression in Eq. (|2.60|) becomes 

-p 2 1 
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(z, -) = -FC A z{\ - z) [2 + e[ln(l - z) + In z] 



For the contribution due to the NGI operator we get (see Eq. ( |2.61|) ) 
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(B17) 



Finally the contribution due to the breakdown of the Ward identity is (see Eq. (|2.62p ) 
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